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leading-order calculations of the QCD background. We apply appropriate cuts to satisfy 
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1. Introduction 

Models with extra dimensions have been proposed to address the hierarchy problem. The 
large extra dimensions model introduced in |l], |2[ suggests that the Standard Model (SM) 
particles live in the usual 3+1-dimensional space, while gravity propagates in a higher 
D-dimensional space. The weakness of gravity with respect to other forces is imposed by 
the large size of the compactified extra dimensions. In this model Newton's constant is 
expressed as 

G^i = Ml = 87ri?"M2+", (1.1) 

where Mp is the Planck mass, Md ~ TeV is the fundamental mass scale, n is the number 
of extra dimensions and R the radius of the compactified space, assumed to be a torus in 
this model. The size of the extra dimensions is determined by their number n, and the 



fundamental mass scale M^, through Eq. |l.l| . For a fundamental scale of the order of 



1 TeV and one extra dimension, its size would be of the order of the solar system size. As 
this would cause large deviations from the inverse square law of gravitation, one simple 
extra dimension is already ruled out by experiments. For more extra dimensions, the size 
gets rapidly smaller, not contradicting submillimetric gravity measurements |^]. 

In the ADD model Q , the graviton corresponds to the excitations of the D-dimensional 
metric. These can be expressed as a tower of Kaluza-Klein (KK) modes. The interaction 
Lagrangian for gravitons and the SM fields is given by 
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with Mp = Mp/v Svr the reduced four-dimensional Planck mass, T^^/ the energy momen- 
tum tensor of the SM fields and n = (ni, n2, n„), a n-dimensional vector of integers, 
labelling the massive gravitons. The interactions of these massive KK states with the 
Standard Model particles are summarised in a set of Feynman rules presented in @, |5[ . It 
is noted that this is an effective linearized theory of interactions valid in the cisplanckian 
region, that is when the centre-of-mass energy of the parton-parton collision is smaller than 
the fundamental scale Md- 

The KK resonances have masses equal to m(^) = \n\/R. This results in a small mass 
gap of order R~^, e.g. for one extra dimension of size 1 //m the mass gap is C(l eV). This 
renders different masses practically indistinguishable and permits replacing the sum over 
discrete mass values by an integral over a continuum with a given density of states. This 
density of states is obtained by considering the number of modes with KK index between 
\n\ and |n| + dn, given by 

dN = Sn-i\n\"-~^dn = Sn-i—oTzrm^^^dm with Sn-i = ^, (1-3) 

M^"- r(n/2) 

Thus the density of states is given by 

r . dN Up 

p^^^ = d^='--'l^'^ ■ ^'-'^ 

As seen from the Lagrangian in Eq. graviton interactions are suppressed by inverse 
powers of Mp. Nevertheless, the summation over the large number of accessible KK modes, 



that is the integration over the density in Eq. 1.4, cancels the dependence on Mp. As an 



example, the inclusive graviton production cross section am is expected from the graviton 
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couplings to be proportional to Mp , which combined with the density in Eq. 1.4 exactly 
cancels the dependence on Mp. This leads to an effective interaction suppressed by inverse 
powers of the fundamental mass scale Mp), thus giving observable effects for Mp near the 
TeV scale. The result is expected from the higher dimensional field theory perspective. 

Interesting phenomenological implications of this extra dimensions model have been 
studied in the literature over the last decade. The processes studied include real graviton 
emission and virtual graviton exchange. A set of processes, such as graviton plus gauge 
boson production, was studied in [Q, |5|, where the Feynman rules were first presented. The 
experimental signature for graviton production is missing energy, as decay into SM particles 
is suppressed by a factor of 1/Mp, which is not compensated for by phase space. Therefore 
gravitons behave as heavy and stable particles, once produced. This missing-energy signal 
does not correspond to a fixed invisible-particle mass as the graviton has a continuous mass 
distribution. This differentiates graviton searches from other Beyond the Standard Model 
(BSM) physics searches, such as supersymmetry which would also give a missing energy 
signal. The experimental signatures depend strongly on the number of extra dimensions, 
n and the fundamental mass scale, Mp, so in principle one should be able to determine or 
at least constrain both of these model parameters. Virtual graviton exchange effects will 
also cause deviations from SM predictions for fermion and boson pair production. Single 
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graviton exchange is an observable not fully calculable. The cross sections are divergent at 
tree level, and thus one has to introduce an ultra-violet (UV) cut-off that is usually taken 
to be of the order of the fundamental scale The search for virtual graviton effects is 
complementary to the production search which is reliably calculated when E <^ M^, and 
will shed light on this cut-off and its dependence on the fundamental scale and the number 
of extra dimensions. 

Both the LEP and Tevatron experiments have set limits on Md. For a brief summary 
of the experimental searches and the current constraints on the fundamental scale and the 
number of extra dimensions see and references therein. A more recent study of the 
dijet angular distribution by the DO experiment at the Tevatron was performed in 0. The 
most stringent current constraints on the fundamental scale are set at 1 — 1.5 TeV but 
these also depend on the number of extra dimensions. The LHC is expected to probe extra 
dimensions effects in all channels mentioned above up to higher scales, due to the larger 
centre-of-mass energy. Recently, a first study of the bounds set on the cut-off scale from 
the dijet cross section measured at the LHC was presented in which are in the few TeV 
region. 

The predictions of this effective theory of graviton interactions fail as we approach 
the quantum gravity scale. As we reach the planckian region where ^/s ~ Md, we need 
theoretical input from quantum gravity, and so it is presently impossible to reliably predict 
experimental signals. An important feature of the present hadron colliders is that the range 
of accessible energies potentially allows the simultaneous probe of cisplanckian, planckian 
and transplanckian regions. Eagerly awaited experimental results from the LHC and future 
colliders could thus provide crucial information on the quantum gravity region. 

The transplanckian region, y/s ^ M^, can also be studied in a fairly model indepen- 
dent way |, 0. In this high energy limit, the two experimental signatures are fermion 
pair production from elastic parton scattering and black hole production |l2|. Calcula- 
tion of the characteristic scales of the system for elastic cross section in the transplanckian 
region [^] shows that it can be described by classical physics. Furthermore, short distance 
gravity effects such as black hole formation are suppressed by considering scattering with 
large impact parameters. For large impact parameters the curvature is small and so we are 
in the limit of weak gravitational field. The scattering amplitude can be computed using 
the eikonal approximation, valid when the scattering angle is small. In this approach, the 
divergent nature of the theory is overcome by resumming an infinite number of ladder and 
cross-ladder Feynman diagrams which give the leading contribution to forward scattering 
1 14]. This has been studied in jl^. Black hole production is also briefly studied in |13| 
in this kinematical region, within a dimensional argument approach. High energy gravita- 
tional scattering has also been studied in other papers in an effort to provide information 
to fill the phase-diagram of different high-energy scattering regimes. In this energy and 
impact parameter space different descriptions such as eikonal scattering. Born scattering 
and string excitations apply in different regions. For a discussion see and references 
therein. 

In this paper, we focus on the transplanckian scattering studied in |13]. In Section 2 
we summarise the calculation of |13| for transplanckian scattering at the LHC, presenting 
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the final result for the eikonal amplitude. We then study the signal and background for 
this process at LHC energies in Section 3. In Section 4, we investigate the extent of the 
computability region by considering the corrections to the eikonal amplitude. In Section 
5, we consider the effect of a possible LHC upgrade to 33 TeV centre-of-mass energy on 
the signal to background ratio as well as the calculability region and compare with 14 TeV 
before we conclude. 



2. Scattering amplitude in the eikonal approximation 

Here we consider graviton exchange in parton-parton scattering in the transplanckian limit. 
The scattering at transplanckian energies can be described by classical physics in a non- 
perturbative approach. In this limit, the gravitational field is weak and non-linear graviton 
couplings can be ignored. In order to tackle the non-perturbative nature of the problem, 
the eikonal approximation is employed. The approximation is valid when the momentum 
transfer is small, or equivalently the scattering angle is small. An infinite set of Feynman 
diagrams is then resummed, and the result is valid when 

< 1 and - t/s < 1. (2.1) 

Details of the calculation are given in |13], where it was found that the eikonal ampli- 
tude is 

Aeik = -2is J (fbe''^^\e''' - 1), (2.2) 
with the eikonal phase x given by 



and 

■(47r)t-isr(n/2)T^/" 



2M]^+2 



(2.4) 



This final result for the amplitude of the infinite series of ladder and cross-ladder diagrams 
can be written more simply as 

A = AnsblFnibcq), (2.5) 
where q = \q±\ ~ -v/— i, and the functions -F„ are 

/•OD 

Fn{y) = -i dxx,h{xy){e'^-'^ -I), (2.6) 







with Jo the zero-th order Bessel function. These functions are calculated numerically to 
high precision in our analysis. 
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3. Phenomenology of transplanckian scattering 



3.1 Signal and background at LHC energies 



The amplitude in Eq. |2.5| was implemented in a Monte Carlo integration program to cal- 
culate the total cross section for this scattering process. For the LHC one should consider 
all possible combinations of quarks and gluons. The matrix element is the same for all 
combinations as expected from this classical picture where the spin of the scattered parti- 
cles does not affect the interaction. The kinematical regions of interest within the eikonal 
approximation are those of \/I/Mo ^> 1 and —i/s <^ 1. This implies that the appropriate 
regions are those of high rapidity separation between the two jets, At/ ^ 1, where 



Ar? = ys - 2/4 = ln(-s/t - 1) 



In 



1 + cos6' 
1 — cos^ 



(3.1) 



with 9 the centre-of-mass scattering angle, and the region of high centre-of-mass energy 
Mjj = VI:$> Md. 

The differential cross section for partonic scattering is 



da nbtse^'' 



dAi] (1 + e^'iy 



\Fniy)\' 



(3.2) 



This must be convoluted with the parton distribution functions of the proton to give the 
hadronic cross section: 



^ dxidx2y^Ji{xi,Q^)fj{x2,Q^) 



dAr] 



dAr] 



(3.3) 



The factorisation scale Q for the parton distribution functions (PDFs) is chosen following 
the prescription in [^], that is = b~'^ li q > b^^ and = otherwise, with bg = 
bcin/qbc)^/^''+^\ The cross sections are calculated with the leading order (LO) PDF set 
MSTW2008LO |16]^. To ensure we stay in the region of validity of the approximation given 
above, suitable cuts must be applied to the final-state jets. The cuts used are px >100 GeV, 
|?7| < 5 for detector acceptance and Mjj > 9 TeV. We note that the detector acceptance 
cuts pt >100 GeV and jr/l < 5 are applied in all subsequent plots in this paper. 

The background for this scattering process is the dijet cross section from QCD scat- 
tering of quarks and gluons. The LO matrix elements squared can be found in The 
background cross section is also calculated with our VEGAS Monte Carlo integrator with 
factorisation and renormalisation scales set equal to pt/2. We note that interference effects 
between the LO QCD background and the signal are negligible, as they are only relevant 
for a small subset of QCD s-channel scattering diagrams. Interference effects would also 
come into play for more combinations of partons if we consider the NLO QCD background. 

The variable chosen for the plots is |A?7|, as experimentally we cannot distinguish be- 
tween the two jets. Therefore events with large momentum transfer, where the partons are 



^We note that in the PDF set used is the leading order CTEQ5L set. We checked that we agree 
by using CTEQ5L to reproduce some of the plots in [|l3[. 
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Figure 1: Differential cross section for the signal for n =6, Md—1.5 and 3 TeV and the QCD 
background. 



scattered backwards retracing their paths, cannot be distinguished from small momentum 
transfer events. Small and negative Ar] regions fall outside the validity of the eikonal ap- 
proximation but their contribution is expected and checked to be small within the eikonal 
formalism due to the small values of Fn (y) at large y = bcQ and thus large momentum trans- 
fer. The differential cross sections of the signal for n = 6, for = 1.5 TeV and 3 TeV and 
the background are shown in Fig. || for ^/s = 14 TeV at the LHC. For Mo = 1.5 TeV, we 
decompose the signal differential cross section da/d\Ar]\ into contributions from positive 
and negative values of Ar/ in Fig. |2[ which confirms that the contribution included from 
regions that fall outside the eikonal approximation is small. 

Several features emerge from Fig. ||. The fall-off at high rapidity separations exhibited 
by both the background and the signal is caused by the decrease of the PDF values at high 
X as the imposed pT cut forces a large dijet mass which implies large momenta fractions 
X and correspondingly small PDFs. For most values of Ar; the values of xi,X2 remain 
close to \/ Smin/s, but at high A77, xi and X2 rapidly approach one. In any case, there is a 
maximum value of allowed Ar/ from kinematical restrictions. The momenta fractions are 
given in terms of y = (ya -|- y4)/2 and Ar] by 

_ _P2LggpQg]^__^ g^j^^ ^2 _ _PZle~^cosh-^. (3.4) 
Multiplying the two, we obtain 

X1X2 = ^7-cosh —— , (3.5) 
2 

which gives a maximum value of Ar/=9.88 for ^/s=14: TeV at the LHC and minimum jet 
transverse momentum pT = 100 GeV. We checked that at 14 TeV the imposed cut on the 
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Figure 2: Contributions to the differential cross section for n =6 and Md=1.5 TeV from positive 
and negative A77. 
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Figure 3: Fraction of total cross section contributed by quark scattering for the signal for n = 6 
at A/u = 1.5 and 3 TeV and the background. 



individual rapidities of the jets |r/| <5 does not alter the behaviour of the cross section in the 
region of high rapidity separation. In this region of high rapidity separation between the 
jets, the QCD background receives significant contributions from BFKL dynamics [19, 20 1, 
which are not taken into account here. 



-7- 



We also note that for this high Mjj cut, the values of xi and X2 are high and the 
dominant contribution to the cross section comes from quark-quark scattering. This is 
shown in Fig. |^ where we plot the fraction of the total cross section we obtain when we 
turn off the gluon contributions. The quark contribution for the signal is larger as all 
combinations of partons have equal weight within the eikonal approximation, while for the 
QCD cross section different colour factors change the weight of each contribution, i.e. the 
effective PDF is g + | E(9 + ^) for the QCD cross section ||T^, compared to 5 + + g) 
for the signal. The values of xi and X2 probed are roughly constant in the range |A?7| < 9. 
The difference between the two values of is caused by the dependence of the PDF 
scale. For the QCD background in Fig. ^, the variation of the quark scattering fraction 
comes from the \^r]\ dependence of the renormalisation and factorisation scales. 

Another effect evident in Fig. |l] is that the differential cross section is larger for smaller 



Mf). From Eqs. 2.4 and 3.2, we know that the differential cross section is proportional to 



^-4(l+2/n)^ The cross section also depends on Md through the y = hcQ argument of the 
functions Fn{y). Large Md and thus small be tend to increase the value of F„ and thus of 
the cross section, but for n = 6 this has a smaller effect than the M^^^^^ factor from 6^ in 
the expression for the cross section. 

The peak structure observed in the differential cross section is due to the form of the 
functions Fn{y) and corresponds to a diffraction pattern of the scattered particles. From 

- (3.6) 



dAT] (1 + e^iy 
the maxima in the partonic cross section are obtained from 

\Fn[y)\ dy 

as y = bcV§/Vl + e'^^. This cannot be solved exactly but the solutions can be estimated 
numerically to locate the peaks, if those are at large A77, so that one can ignore the second 



term in the right-hand side of Eq. |3.7| . The number of extra dimension changes the structure 
of the functions Fn- For comparison, in Fig. ^ we show the differential cross section for 
n =2, 4, 6 on a linear scale at 14 TeV with the same cuts as in Fig. ||. The different position 
of the peaks for each n could, at least in principle, be used to experimentally extract the 
number of extra dimensions. 

3.2 Comparison of eikonal signal to other signals 

In order to better understand the shape of the cross section, we consider a test case where 
FQ^y) is fixed at its y = value, i.e. its value at large rapidity separations. In this case, the 
cross section is expected to keep increasing monotonically on moving from large to small 
rapidity separations, giving a straight line at sufficiently large Arj in a log plot. This is 



explained simply by Eq. 3.2, which at large Ar] simplifies to 

da 



7rbtse-^^\Fn{y)\'. (3.8) 
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Figure 4: Differential cross section for n =2, 4, 6 and Mjj—l.b TeV. 

Note that similar behaviour is expected for a four-particle contact interaction such as the 
one proposed in p^], where the contact interaction for left-handed quarks has the form 



nog 



= Yf^^L'y^'lLqLltJ.qL with no = ±1. 



(3.9) 



The amplitudes squared are given in [^Tj. For comparison, we show in Fig. I the cross- 
section results for the signal, for fixed F^d/) = Fq{0), for the contact interaction and for the 
functions e"^'' and e^^/(l + e^^)^ to illustrate the different functional behaviours. The 
contact interaction constant j IS? and the coefficients of the two functions are adjusted 
so that their values match those of the signal at the matching region of large 

Of course, there are strong experimental limits on the contact interaction strength from 
several experiments, see for example the summary table in [p2| and the recent results from 
the DO collaboration 0]. Recently, the first LHC results were presented by the ATLAS [23| 
and CMS collaborations. The two collaborations searched for deviations from the SM 
prediction for the angular distribution of dijet events and raised the 95% c.l. exclusion limit 
to 3.4 and 4 TeV respectively. 

Here, we are only interested in establishing the functional dependence of the differential 
cross section on A// for the different cases listed above. We note that the behaviour 
described by Eqs. p.2| , p.8| refers to the partonic cross section, but it is also shown by the 
hadronic cross section for all but very high Ar/. This can again be explained by the high 
Mjj cut and the fact that the cross section is a rapidly decreasing function of Mjj. This 
leads to approximately constant values of x\ and X2 and thus s, which do not allow the 
integration in Eq. |3.3| to alter the functional dependence of the cross section on Ar/. 



-9- 



100 



42 



10 r 



1 r 



< 0.1 

b 



0.01 r 



0.001 r 



0.0001 



/5 = 14 TeV, Mjj > 9 TeV 

r^. — I 1 1 1 1 1 1 r~ 

- Signal A/d=1.5 TeV 

F(y)=F(0) 

Contact Interaction 

--l^V, (l+exp(A,7))2 
■ - . ^ 




_L 



_L 



4 5 
\Arj\ 



Figure 5: Differential cross section for different cfioices of interaction. The contact interaction 
constant has been adjusted to be of the same order of magnitude in the matching behaviour region, 
as are the coefficients of the two functions. 



We also note that the differential cross section is symmetric in Arj for the contact 
interaction and for the case of fixed F(y), which leads to a factor of two difference at the 
matching region between the full treatment and this simplified test case of constant F. 
The graviton signal cross section is not symmetric in Ary, as the value of y increases with 
decreasing Ar] and Fn{y) — )• as y — )• oo. On the other hand, the QCD background is 
symmetric in Arj, assuming a renormalisation and factorisation scale choice symmetric in 
Arj. 



In addition to the dimension-6 operators of Eq. |3.9| , which arise in BSM theories such as 
technicolour, but also in extra dimensional theories from one- loop graviton exchange |25|, 
one can also compare the transplanckian signal to the dimension-8 operators that arise in 
extra-dimensional theories. In this case, the amplitude for a single virtual graviton ex- 
change is represented by the energy momentum tensor squared and the UV divergence is 
regularised arbitrarily by introducing the coupling constant A^^. This is used to calcu- 
late the virtual graviton effects in the cisplanckian region. The value of the constant is 
expected to be nearly the value of Md but precise numerics are unknowable. The corre- 
sponding matrix elements for dijet production have been calculated in [26|. The differential 
cross section da/dAr] for the signal is added to the background for the dimension-6 and 
8 operators and shown in Fig. ^ for Aj'=3 TeV. We note that we choose the value of the 
fundamental scale M^, the strength of the contact interaction A and the UV cut-off At to 
be the same for illustration purposes. One can appropriately rescale the cross sections for 
the dimension-6 and dimension-8 operators, as in the region where the signal dominates 
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Figure 6: Differential cross section for interactions originating from dimension-8 operators in com- 
parison with dimension-6 operators, the transplanckian signal for n = 6 and the QCD background. 

over the background, they have a simple dependence on the coupling constant: A""^ and 
A^*^ respectively. 

The shape of the signal can evidently be used to distinguish between the different 
types of interactions with the eikonal amplitude giving a distinctive peaking structure. We 
emphasize here that the calculation for the dimension-6 and dimension-8 operators signals 
is strictly not valid in this region of high centre-of-mass energy larger than the cut-off scales. 
The purpose of Fig. ^ is simply to show the functional behaviour of the cross section for 
the chosen region of Ar] for different operators. 

3.3 PDF and PDF scale dependence of the cross sections 

The dependence of the results on different renormalisation and factorisation scale choices 
has also been investigated. The scale affects the signal through the PDF values and the 
background through both the PDF values and the explicit value of Og in the matrix el- 
ements, as we have set the renormalisation scale equal to the factorisation scale in this 
study. The relations connecting Ar/ and possible scale options are 



PT 



1 + e^'' ' 



n 

qbc 



-l/(n+l) 



(3.10) 
(3.11) 
(3.12) 



with be given in Eq. 2.4. 
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Figure 7: Dependence of the cross section on scale choices for both the background and the signal 
for n — 6. 



The effect of the scale choice on the LO differential cross sections is shown in Fig. |^. 
The choice Q = q = \f—i matches the prescription in [^] for At/ larger than ~7, in 
agreement with the crossing between q and 1/fec- In general, we see that smaller scale 
choices lead to larger cross sections. One should of course consider next-to-leading order 
calculations to reduce this dependence on the scale choice, as has been pursued in [^]. 

As an additional check, we also consider two different PDF sets. The differential cross 
section for the signal at Mjj = 3 TeV and the background is shown for two different leading 



order PDF choices, MSTW2008LO [|T§ and CTEQ6L |2|], in Fig. |. There is a difference 
of about 35% between the two sets for the signal, while the difference is smaller for the 
background. For the signal the difference is mostly caused by the difference of the PDF 
values for the valence u and d quarks which at x = 0.7 differ in the two sets by almost 
20% and 40% respectively at Q = 100 GeV. For the background, this difference is partially 
compensated to some extent by the value of , for which the values of the two sets differ 
in the opposite way by about 20% at Mz- In Fig. we show the signal to background 
ratio for the two PDF sets. The difference in the signal is partially compensated by the 
difference in the background for small |A?7| values. 

3.4 Total cross section for eikonal signal 

In Fig. 0, the total cross section is shown as a function of the Mjj cut for both the signal 
and the background. The rapidity separation is forced to be in the range 3 < |At7| < 4, 
and the px > 100 GeV and \rj\ < 5 cuts remain in place. The choice of the acceptance 
region of \Ar]\ is a compromise between staying in the kinematical region allowed within the 
eikonal approximation (large Arj) and maximising the signal to background ratio (small 
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Figure 8: Dependence of the cross section on 
PDF clioices for both the n = 6 signal and the 
background. 



Figure 9: Dependence of the signal to back- 
ground ratio on the choice of PDF set, for the 
signal of ri = 6. 



Ar/). We note that the background falls more rapidly than the signal with increasing 
Mjj". Therefore the search for signal against background will be easier for a higher cut 
provided the event rate remains high enough for detection. We also note that the cross 
section is smaller for larger M^), as it is suppressed by inverse powers of M^. Moreover 
the cut needed to ensure we remain in the transplanckian region increases with larger Mq, 
leading to smaller cross sections and lower event rates. We thus expect better prospects of 
discovery for smaller M£, values. 



Finally, in Figs. |11| and we show the dependence of the total cross section on 
the fundamental scale M/j, fixing the Mjj cut at 2>M£, and the more conservative QMd 
respectively, again for 3 < |Ar/| < 4. Given the current experimental limits on the fun- 
damental scale Md with the centre-of-mass energy of 14 TeV available at the LHC, we 
can barely reach the transplanckian region for the more conservative Mjj=QM£) choice, 
with the PDF values suppressing the cross sections. Although the more reliably calculated 
limits from graviton emission are still well under 2 TeV, the indications from dimension-8 



and dimension-6 induced operators, such as those discussed in [25| are that Mo may need 
to be higher, of order several TeV. Although such estimates are uncertain and dependant 
on the yet-unknown underlying theory of quantum gravity, if correct they would put the 
transplanckian region nearly out of reach for 14 TeV LHC. However, the recently discussed 
higher energy LHC of 33 TeV could still reach transplanckian validity in that case. That 
will be the focus of the penultimate section, but first we discuss how to quantify roughly 
the calculability of the transplanckian signal in the next section. 

4. Computability region of transplanckian scattering signal 



The conditions for the calculation to be valid are given in Eq. 2T. These relate to the 
constraint on the energy to be in the transplanckian region and the small scattering angle 
for the eikonal approximation to be valid. In the transplanckian region scattering can be 
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Figure 10: Total dijet cross section as a function of tlie Mjj cut for tlie signal at Mjj —1.5 and 
3 TeV, for n — 6 and the background. 
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Figure 11: Total cross section for n = 6, as a Figure 12: Total cross section for n = 6, as a 



function of Md for A/™" = 3Md 



function of Md for M™" 



6Md. 



described by classical physics and is thus only calculable by non-perturbative methods. 
Imposing also the constraint of large impact parameter, so that the gravitational field is 
weak, allows us to use the eikonal approximation to calculate the scattering amplitude. The 
eikonal approximation is only valid for small scattering angles which imposes the second 
condition. Corrections to the eikonal amplitude originate from both classical and quantum- 
mechanical effects. The classical effects originate from neglecting terms of order —t/s and 
Rs/b which are assumed to be very small. The emergence of the terms is discussed below. 
For a collision with centre-of-mass energy ^/s, the Schwarzschild radius of the system 
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is formed as: 



Rs = ^ 



{GdV~s)'/'^'''-'\ (4.1) 



(n + 2) 

with Go the D-dimensional Newton's constant. In the region where ^/s ^ q ^ b^^, b < be 
and the eikonal phase is large, and we are in the classical region. By dimensional analysis 
the scattering angle of a collision is related to the impact parameter b and the Schwarzschild 
radius by ~ {Rs/b)^^^ . Therefore the expansion parameter —t/s is related to the impact 
parameter of the collision (ignoring factors of order one) by 



Therefore in this region the two conditions coincide. A small —t/s also ensures that the 
impact parameter is larger than the Schwarzschild radius so that the gravitational field 
is weak and the risk of black hole formation is avoided. By including non-linear effects 
Eq. |4.2| will be modified by a factor oi \ + 0{t/ s). 

The other aspect that needs to be controlled for this semi-classical calculation to be 
valid is the quantum nature of the particles. When q <^ b~^ , b ^ b^ the eikonal phase is of 
the order one and quantum effects are important. In this region the size of the correction 



for small t is [13|: 



1 fRsV'-'-' [Ml^ 



sbl \ be J \ s 



Combining the two, the leading classical gravity corrections to the eikonal amplitude are 

(4.4) 



o{ --\+o 





t 




e = 








s 





Ml 
s 

In addition to these, there are quantum gravity corrections to the eikonal amplitude. These 
may be model-dependent string corrections which we do not consider here. 

In order to visualise and quantify the computability of the two-to-two scattering, we 
introduce the parameter e, defined as 

rf) ■ C'^) 

For any assumed or given value of the parameter e is well determined from kinematical 
quantities we used in the previous sections to study the two-to-two scattering, as there is 
direct correspondence between —t/s and the rapidity separation between the two jets Ar/ 
through 

-- = TT^- (4.6) 

s 1 + e^^ 

Moreover, -y/i is given by the dijet mass. The dependence of e on s and Ar] for n = 6 is 
shown in the contour plot of Fig. [l^. The lines are contours of constant e with the value for 
each colour line given in the legend. From this contour plot we can determine the region 
in the Mjj and Ar] space where the calculation is reliable. Computability requires a small 
value of e which corresponds to large Ar; and large Mjj/Mjj- In the following section, we 
will use e to visualise the extent of the calculability region of the theory. 
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Figure 13: Dependence of e on Arj and Mjj/Md for n —6. The lines are contours of constant e. 
5. Comparison of HE LHC and LHC 

Increasing the centre-of-mass energy of the colUder would enable us to probe collisions with 
partonic centre-of-mass energy much greater than Md and thus well into the transplanckian 
regime. The proposed energy upgrade for the LHC (HE LHC) based on 20 T magnets, 
would increase the beam energy to 16.5 TeV |^^. For the HE LHC, the higher centre-of- 
mass energy of 33 TeV allows us to set the Mjj cut higher while at the same time the 
total rate for the signal is increased significantly. For comparison with Fig. ^ we show the 
corresponding plot for ^/s= 33 TeV in Fig. |l^ We note that in this plot the cross section 
units are nb compared to pb we had in Fig. 0. Adding all cases in one plot to better 



visualise the different effects present, we show in Fig. 15 the differential cross section for 
M£) = 3 TeV and n=2, 4, 6 with the same Mjj cut. While the signal to background ratio 
is similar for both centre-of-mass energies, the cross sections are four orders of magnitude 
larger for =33 TeV leading to more easily measurable event rates. 

Moreover, we now have the chance to increase the Mjj cut retaining measurable cross 
sections. The effect of increasing the dijet mass cut is shown in Fig. |l^. On increasing 
the dijet mass cut the crossing point of signal and background moves to the right, and 
therefore the region where the signal is larger than the background is extended to higher 
Ar] values. The importance of this is explored in the rest of the section where we examine 
the effect of increasing the centre-of-mass energy on the calculability region. 

On increasing the centre-of-mass energy, the calculability region gets wider as shown by 
the value of the parameter e as a function of kinematics. Here, we calculate the differential 
cross section da/dArj as a function of Ar] using the parameter e to visualise where the 
amplitude is reliably calculable. The pT and r] cuts of Section 2 are imposed here to account 
for detector acceptance. We note that in this section we restrict ourselves to calculating 
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Figure 14: Differential cross section da/d\lS.r]\ at =33 TeV for =1.5 TeV for different n. 

Md=?> TeV, Mjj >9 TeV 




|Ar?| 

Figure 15: Differential cross section da / d\/S.i-]\ at yfs =14 and 33 TeV for Md =3 TeV for different 
n. 

the signal and background for positive values of Ar;, rejecting the region of Ary <0 which 
corresponds to large momentum transfer scattering which cannot be calculated within the 
eikonal approximation. The contribution of this region is expected to be sufficiently small. 
The value of e is calculated at each value of A?/ using the invariant dijet mass cut squared 
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Figure 16: Differential cross section da/d\Ari\ at -^i =33 TeV for Mjj —3 TeV and n=6 varying 
the Mjj cut. 



in place of s in Eq. calculated for Md=1.5, 3, 5 TeV with M^j" = 9 TeV at the LHC 
in Fig. |l^. The signal lines have three colours depending on the value of e: green (solid) for 
e <0.15, blue (dashed) for 0.15 < e < 0.3 and red (dotted) for 0.3 < e < 0.5. The plots do 
not extend into regions where e > 0.5 as there the calculation is not reliable. Calculability 
improves (i.e., e becomes smaller) by increasing the centre-of-mass energy as we can then 
probe deeply into the transplanckian region. This is illustrated in Fig. where the dijet 
invariant mass cut has been raised to 15 TeV, an energy not accessible at present LHC 
energies. Now the M£)=5 TeV line has turned green (e < 0.15) which shows that higher 
Md scales can now be probed reliably, with the signal to background ratio being larger 
and the green (solid) region extended for all M^. 

In a similar way, computability can be visualised for the differential cross section 
da/dMjj. Keeping the same px and individual jet ?? cuts of Section 3, we calculate 
da/dMjj for two different At] bins: 2 < At] < 4 and 4 < At] < 6. Here e is calculated using 
the minimum value of Arj in each bin to give the most stringent limit on computability. 
The cross section is calculated for M/)=1.5, 3, 5 TeV at both 14 and 33 TeV for the same 



range of Mjj values and shown in Fig. 19 and Fig. respectively. We note that using 
the minimum value of A77 in each bin to calculate e for da/dMjj gives generally stronger 
limits for computability compared to using the Mjj cut for da/dArj, as da/dMjj is a 
rapidly decreasing function of Mjj whereas da /d Arj is not a decreasing function of Arj. 
In all plots we show the QCD background, as we need to keep in mind that in addition to 
computability, the signal to background ratio is crucial for any experimental search. In the 
plots, we see again that setting a higher Arj cut improves the reliability of the calculation 
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Figure 17: Differential cross section da/dAii for Figure 18: Differential cross section da/dArj for 
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Figure 19: Differential cross section da/dMjj for n = 6, as a function of Mjj for 2 < A77 < 4 
and 4 < A?7 < 6 at TeV. 



but decreases the signal to background ratio. Thus, an intermediate region of At/ is ideal. 
Again increasing the centre-of-mass energy gives better signal to background ratio and a 
wider calculability region. 

6. Conclusions 

We have studied transplanckian scattering at the LHC within the eikonal approximation. 
In the kinematical limits of high rapidity separation between the two jets and high centre- 
of-mass energy, we expect the distinctive features of the signal to allow us to distinguish it 
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Figure 20: Differential cross section da/dMjj for n = 6, as a function of AIjj for 2 < Arj < 4 
and 4 < Ar; < 6 at Vi=33 TeV. 



from the QCD background, which we have also calculated and presented. At the 14 TeV 
LHC, given the current limits on the fundamental mass scale, it is only marginally possible 
to reach the transplanckian regime. We have studied how the signal and background are 
affected by the fundamental mass scale, the scale choice for the PDFs and the approach 
needed to maximise the signal to background ratio. We also considered the computability 
of the eikonal signal and considered the corrections to the eikonal amplitude to establish 
the region of parameter space where the eikonal amplitude is most reliably calculable. 

We considered how a possible LHC energy upgrade to 33 TeV centre-of-mass energy 
will improve the prospect of observing transplanckian scattering at the LHC. Our studies 
have shown that even with > 5 TeV there is in principle a calculable transplanck- 
ian kinematic region with good signal rate (~ fb) and favourable signal to background 
ratio. This constitutes another illustration of how increasing the energy is typically more 
important in most "gravity phenomenology" than increasing luminosity. Thus, from this 
perspective the HE LHC presents qualitative improvement in covering the parameter space 
of extra dimensional theories with its prospects for probing the cisplanckian, planckian and 
transplanckian regimes, whose signals could be determined through control of kinematic 
variables. 
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